Searching for dark matter sterile neutrino in laboratory 
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If the dark matter of the Universe is made of sterile neutrinos with the mass in keV region they can be 
searched for with the help of X-ray satellites. We discuss the prospects of laboratory experiments that can be 
competitive and complimentary to Space missions. We argue that the detailed study of j3 decays of tritium and 
other nuclei with the help of Cold Target Recoil Ion Momentum Spectroscopy (COLTRIMS) can potentially 
enter into interesting parameter range and even supersede the current astronomical bounds on the properties of 
dark matter sterile neutrino. 
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I. INTRODUCTION. 

The nature of Dark Matter (DM) in the Universe is a puzzle. 
Many different hypothetical particles coming from physics 
beyond the Standard Model (SM) were proposed to play a role 
of dark matter particle; none of them have been discovered 
yet. In this paper we will discuss possibilities for a laboratory 
search of one of the dark matter candidates: sterile neutrino 
with the mass in keV region. 

In short, here is the case for sterile neutrino as a dark mat- 
ter particle. There are not that many experimental facts in 
particle physics which cannot be described by the Standard 
Model. These are neutrino oscillations (neutrinos of the SM 
are exactly massless and do not oscillate), dark matter (the SM 
does not have any stable neutral massive particle) and baryon 
asymmetry of the Universe (substantial deviations from ther- 
mal equilibrium, needed for baryogenesis, are absent for ex- 
perimentally allowed mass of the Higgs boson; in addition, 
it is a challenge to use CP violation in Cabibbo-Kobayashi- 
Maskawa mixing of quarks to produce baryon asymmetry in 
the SM). This calls for an extension of the SM. Perhaps, the 
most economic one, that can describe all these phenomena in 
a unified way is the vMSM of [Ij, i2i]. In this model three 
leptonic singlets (other names for them are right-handed, Ma- 
jorana or sterile neutrinos) are added, making the structure of 
quark and lepton sectors of the theory similar, up to introduc- 
tion of Majorana mass for the new leptonic states. The Ma- 
jorana nature of singlet fermions leads to non-zero masses for 
active neutrinos and, therefore, to neutrino oscillations, solv- 
ing in this way one of the SM problems. The lightest of these 
new particles with the mass in keV region can have a life- 
time greater than that of the Universe ^ and thus can play a 
role of (warm) dark matter [4]. The preference for keV mass 
scale is coming from the cosmological structure formation ar- 
guments related to the missing satellites problem jsl 0] and 
to cuspy DM distributions in cold dark matter cosmologies 
||7|,i8i]. For other astrophysical applications of keV sterile neu- 
trinos see ll9[ [lol[TTl[l2[ll3[l . The presence of two other heav- 
ier fermions with the mass in GeV region leads to gen- 
eration of baryon asymmetry of the Universe |2] via resonant 
sterile neutrino oscillations iH and electroweak sphalerons 
ifisll . These fermions can be searched for in particle physics 



experiments with high intensity proton beams lfT6l[T7ll . 

The only way considered up to now to detect the dark mat- 
ter sterile neutrino is through astrophysical X-ray obser- 
vations mini- The interaction of with intermediate vector 
bosons W and Z and ordinary charged leptons and neutrinos is 
suppressed by the so-called mixing angle Q — mo /Mm, where 
mo and Mm are respectively the Dirac and Majorana masses 
of sterile neutrino. The main (but undetectable) decay modes 
of sterile neutrino are ^ 2v + V, ^ V + 2 v. In addition, 
sterile neutrino has a radiative decay channel A^ ^ V + 7, A^ ^ 
V + 7, producing a narrow line in X-ray spectrum coming from 
dark matter halos of different astronomical objects. Recently 
a number of constraints on the mixing angle Q became avail- 
able, coming from the analysis of X-ray data of Chandra and 
XMM-Newton satellites |19, 20, 2ll|2i|2lf24, 25, 26, 27]. 
It is expected that the best results will come from the analy- 
sis of dwarf satellite galaxies in the Milky way halo, as having 
the largest mass-to-light ratios and weakest X-ray background 

In addition to X-ray constraints there are bounds on the 
mass and momentum of dark matter sterile neutrinos coming 
from analysis of Lyman-a forest clouds and structure forma- 
tion [28, 29j 30]. They depend, however, on specific mecha- 
nism of cosmological production of sterile neutrinos Isill . The 
most conservative limit on the mass, Mn > 0.3 keV is com- 
ing from the anaWsis of rotational curves of dwarf spheroidal 
galaxies (^2. HH [34|] (Tremaine-Gunn bound). 

Yet another constraint comes from the requirement that the 
amount of sterile neutrinos produced in the early universe due 
to the mixing with ordinary neutrinos must be smaller than the 
amount of the dark matter observed. In the absence of entropy 
production due to decays of heavier singlet fermions ll3lll and 
assuming that the standard Big Bang theory is valid at tem- 
peratures below few hundreds MeV these bounds are stronger 
than those coming from X-ray observations for M^v < 3.5 keV 
llssi [36I1 . However, relaxing the above-mentioned assumptions 
can eliminate these constraints I3l, i37ll . 

To summarize, the most conservative constraints on the 
dark matter sterile neutrinos are coming from X-ray observa- 
tions and from rotational curves of dwarf galaxies; only those 
will be used in what follows. In Fig. [T] we present the main 
X-ray bounds taken from ll2ll I25I1 in the mass range allowed 
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by the Tremaine-Gunn bound, for comparison with the pro- 
posed laboratory experiment sensitivity. We stress that these 
constrains are purely observational and do not depend on any 
theoretical bias; the only assumption behind is that the sterile 
neutrinos constitute 100% of the dark matter in the universe. 
If only a fraction p of dark matter is in sterile neutrinos, the 
X-ray constraints are weaker by a factor of p. For p consid- 
erably smaller than 1 the Tremaine-Gunn bound is also not 
applicable. 

Imagine now that some day an unidentified narrow line will 
be found in X-ray observations. Though there are a number 
of tests that could help to distinguish the line coming from 
DM decays from the lines associated with atomic transitions 
in interstellar medium, how can we be sure that the dark mat- 
ter particle is indeed discovered? Clearly, a laboratory ex- 
periment, if possible at all, would play a key role. Current 
bounds in the interesting mass region were mostly based on 
kink search in beta decay, inspired by possible discovery of 
17 keV neutrino. The present bounds |38] are given in Fig.[T] 
and are much weaker than required to compete with X-ray ob- 
servations. Fig. [T] demonstrates that the search for DM sterile 
neutrino in terrestrial experiments is very challenging, as the 
strength of interaction of DM sterile neutrino with the matter 
is roughly times weaker than that of ordinary neutrino. In 
this paper we analyze different processes where DM sterile 
neutrino can be searched for in the laboratory. We argue that 
the only potential possibility is provided by a precise study of 
the kinematics of beta decays. We note, in particular, that this 
kind of study is not impossible in a view of a novel momentum 
space imaging technique (Cold Target Recoil Ion Momentum 
Spectroscopy, COLTRIMS), proposed and developed at the 
end of nineties (for reviews see ll39ll40ll ). 

The paper is organized as follows. In second section we 
will discuss different reactions in which the DM sterile neutri- 
nos can potentially manifest themselves. In third section we 
consider the requirements to the j3 -decay recoil experiments 
that can enter into interesting region of mixing angle for DM 
sterile neutrino and provide arguments that they could be po- 
tentially feasible. The last section is conclusions. 



II. DM STERILE NEUTRINOS IN THE LABORATORY. 

In the vMSM interaction of DM sterile neutrino with 
fermions of the SM can be derived from the standard 4- 
fermion weak interaction by replacements Va ^ Va + daN, 
where a = e, /i, T, 0^=520a, and we work in the lowest or- 
der in mixing angles Qa- Thus, sterile neutrinos participate in 
all reactions the ordinary neutrinos do with a probability sup- 
pressed by 0^. Additionally, the fact that they are Majorana 
particles, — (c is the sign of charge conjugation), leads 
to lepton number non-conservation. 

From very general grounds the possible experiments for the 
search of sterile neutrinos can be divided in three groups: 

(i) Sterile neutrinos are created and subsequently detected 
in the laboratory. The number of events that can be associated 
with sterile neutrinos in this case is suppressed by 0"* in com- 
parison with similar processes with ordinary neutrinos. The 



smallness of the mixing angle, as required by X-ray obser- 
vations, makes this type of experiments hopeless. For exam- 
ple, for sterile neutrino mass nis — 5 keV, the suppression in 
comparison with neutrino reactions is at least of the order of 

io-'9. 

(ii) Sterile neutrinos are created somewhere else in large 
amounts and then detected in the laboratory. The X-ray Space 
experiments are exactly of this type: the number density of 
sterile neutrinos is fixed by the DM mass density, and the lim- 
its on the X-ray flux give directly the limit on 0^ rather than 0^ 
as in the previous case. Another potential possibility is to look 
for sterile neutrinos coming from the Sun. The flux of sterile 
neutrinos from, say, pp reactions is Fat ~ 6 x lO'^S^ /cm^s. 
The only way to distinguish sterile neutrinos from this source 
from electronic neutrinos is the kinematics of the reactions 
VgH pe and Nn — > pe, which looks hopeless. For higher 
energy sources, such as neutrinos, the emission of ster- 
ile would imitate the antineutrinos from the Sun due to the 
reaction Np ne^ which is allowed since is a Majorana 
particle. However, this process is contaminated by irremov- 
able background from atmospheric antineutrinos. Even if all 
other sources of background can be eliminated, an experiment 
like KamLAND would be able to place a limit of the order of 
0^ < 3 X 10^^, which is weaker than the X-ray limit for all 
possible sterile neutrino masses obeying the Tremaine-Gunn 
bound ll32l [33l 13411 . The current KamLAND limit can be 
extracted from IHTIl and reads 0^ < 2.8 x lO^'*. The sterile 
neutrino can also be emitted in supernovae (SNe) explosions 
in amounts that could be potentially much larger than 0^/v, 
where Fy is the total number of active neutrinos coming from 
SNe. The reason is that the sterile neutrinos interact much 
weaker than ordinary v and thus can be emitted from the vol- 
ume of the star rather than from the neutrino-sphere. Using the 
results of |j3|], the flux of SNe sterile neutrinos due to — A^ 
mixing is Fat ~ 5 x 10^0^ {ms/k&W)'^ F^. In spite of this en- 
hancement, we do not see any experimental way to distinguish 
the A^ and (v, v) induced events in the laboratory. 

(iii) The process of sterile neutrinos creation is studied in 
the laboratory. In this case one can distinguish between two 
possibilities. In the first one, we have a reaction which would 
be exactly forbidden if sterile neutrinos are absent. We were 
able to find just one process of this type, namely S invisible, 
where S is any scalar boson. Indeed, in the SM the pro- 
cess 5 — > VV is not allowed due to chirality conservation, and 
5 ^ VV is forbidden by the lepton number conservation. With 
sterile neutrinos, the process S vN may take place. How- 
ever, a simple estimate shows that the branching ratios for 
these modes for available scalar bosons such as k'^ or are 
incredibly small for admitted (by X-ray constraints) mixing 
angles. So, only one option is left out: the detailed study of 
kinematics of different j3 decays. 

An obvious possibility would be the main pion decay mode 
^ jiv with creation of sterile neutrino A^ instead of the 
active one. This is a two body decay, so the energy muon 
spectrum is a line with the kinetic energy (m^j — mp)^ /Im-n = 
4.1 MeV for decay with active neutrino and {{mj^ — m^Y — 
m1)/2m-„; for decay with massive sterile neutrino. Thus, for 
nis of keV order one needs the pion beam with energy spread 
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less than 0.01 eV to distinguish the line for sterile neutrino, 
which seems to be impossible to get with current experimental 
techniques. 

In the case of beta decay there are two distinct possibili- 
ties. One is to analyze the electron spectrum only. In this 
case the admixture of sterile neutrinos leads to the kink in the 
spectrum at the distance nis from the endpoint. However, the 
distinguishing a small kink of the order of 6^ on top of the 
electron spectra is very challenging from the point of view of 
statistically large physical background and nontrivial uncer- 
tainties in electron spectrum calculations. 

The case of full kinematic reconstruction of beta decay of 
radioactive nucleus is more promising and will be analyzed in 
the next section. 



III. COLTRIMS AND /3 -DECAYS. 

The idea of using beta decay for sterile neutrino detection 
is quite simple: measuring the full kinematic information for 
the initial isotope, recoil ion, and electron one can calculate 
the neutrino invariant mass on event by event basis. In an 
ideal setup of exact measurement of all these three momenta 
such an experiment provides a background-free measurement 
where a single registered anomalous event will lead to the 
positive discovery of heavy sterile neutrino. This idea was 
akeady exploited at time of neutrino discovery and testing of 
the Fermi theory of j3 decay ll42ll . It was also proposed to 
use full kinematic reconstruction to verify the evidence for 
17 keV neutrino found in the kink searches [43, 44] to get 
rid of possible systematics deforming the beta spectrum. Re- 
cently, bounds on sterile neutrino mixing were achieved by 
full kinematic reconstruction of ^^'"K isotope confined in a 
magneto-optic trap |45] but for a neutrino in the mass range 
0.7 — 3.5 MeV, what is much heavier than considered here. 
For 370 — 640 keV mass range a similar measurement was 
performed in electron capture decay of ^^Ar l46ll . We will 
discuss below a possible setup for a dedicated experiment for 
a search of keV scale DM sterile neutrino. 

Let us consider an idealized experiment in which a cloud of 
/3 -unstable nuclei, cooled to temperature T, is observed. For 
example, for the normal beta decay is 

^H^^He + e + v^, 

while in presence of sterile neutrino in about 6^ part of the 
events (up to the kinematic factor) the decay proceeds as 

where is a sterile neutrino in mostly right-handed helicity 
state. Suppose that it is possible to register the recoil momen- 
tum of the daughter ion and of the electron with high enough 
accuracy. Indeed, existing COLTRIMS experiments are able 
to measure very small ion recoil 1391 l40ll . They are utilized 
for investigation of the dynamics of ionization transitions in 
atoms and molecules. The ion momenta is determined by time 
of flight measurement. A small electric field is applied to the 
decay region to extract charged ions into the drift region. After 



the drift region the ions are detected by a position sensitive de- 
tector, which allows to determine both the direction of the mo- 
menta and the time of flight. Characteristic energies of recoil 
ion in beta decay is of the order of the recoil momenta mea- 
sured by existing COLTRIMS in ion-atom collisions. Preci- 
sions currently achieved with such apparatus are of the order 
of 0.2 keV for the ion momentum fsl El, gl gl H . 

Electron detection is more difficult, as far as the interest- 
ing energy range is of the order of 10 keV for decay (or 
greater for most other isotopes). This is much higher than typ- 
ical energies obtained in atomic studies. One possible solution 
would be to use the similar time-of-flight technique as for the 
recoil ions, but with adding magnetic field parallel to the ex- 
traction electric field, thus allowing to collect electrons from 
a wider polar angle. In existing applications such a method 
was used for electrons with energies of only 0.1 ke V [51, 52i ■ 
In llssll retarding field was added in the electron drift region 
allowing to work with electrons of up to 0.5 keV energies. Al- 
ternatively, one may try to use electrostatic spectrometers for 
electron energy measurement, as it was proposed in ll43ll44ll . 
On the one hand, the latter method allows to use the electron 
itself to detect the decay moment for recoil time of flight mea- 
surement. On the other hand, it is hard to reach high polar 
angle acceptance with this method, thus losing statistics. 

The decay moment needed for the time-of-flight measure- 
ment can be tagged by registering the Lyman photon emission 
of the excited ion or by the electron detection, if electron en- 
ergy is determined by a dedicated spectrometer. Note, that for 
case Lyman photon is emitted only in about 25% of the 
events 1441 , so the photon trigger also induces some statistics 
loss. 

According to (3^ it is possible to achieve sensitivity for 
measuring normal active neutrino masses of 10 eV for each 
single event; the accuracy needed for the case of sterile neutri- 
nos is considerably less than that as the mass of is expected 
to be in the keV region. Moreover, the measurement in the 
latter case is a relative measurement, which is much simpler 
than absolute measurement of the peak position required for 
active neutrino mass determination. 

Let us estimate the required precision of momentum mea- 
surements and source temperature. Suppose that an initial 
molecule with mass M decays at rest into recoil ion, electron 
and neutrino with momenta p, k and q = p + k respectively. 
The energy release will be denoted by Q. Then the neutrino 
mass can be defined from ion and electron momenta as 

ml^iQ-E,-E„f-{p + k)\ 



where Eg = y/ + k^ — nie and Ep = fvfi + p^ — M are 
electron and recoil ion kinetic energies. It is immediately 
seen, that for measuring keV neutrino mass precision of 
0.5 keV in momenta measurement would be sufficient. For 
example, for decay this means 0.5 x 10^^ precision in mo- 
mentum measurement. 

Let us turn now to the question of the temperature of the 
cloud. Nonzero thermal velocity v of the initial molecule 
spoils the measurement. The measured values of the momenta 
would be p = p + My and k = k + m^v so the measured mass 
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FIG. 1: Constraints on the mixing angle 6 of sterile neutrino with 
active neutrino from X-ray observations of Large Magellanic Clouds 
and Milky Way by XMM-Newton and Milky Way by HEAO-1 satel- 
lites and from kink searches in beta decay. The X-ray bound is in 
the assumption that sterile neutrinos constitute p = 100% of the DM, 
for smaller p the bound is relaxed accordingly. The boundaries of 
the parameter space accessible to j3-decay experiments are shown 
for various values of the source temperature, kinematic cuts and col- 
lected statistics. 
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where q = p + k and we have neglected terms of the order 
of nie, Q, and |k| compared to the ion mass M. The average 
squared thermal velocity is (v 
momenta (in the rest frame) is |q 
that for reasonably low temperatures the last term leads to the 
dominant error. Assuming an isotropic thermal probability 
distribution of decaying atoms P{\) °^ exp(— Mv^/2r) it is 
easy to find (with exponential accuracy) the probability to get 
nonzero value of = nis from an event originating from a 
j3 decay to massless ordinary neutrino, 



P{ms) °= exp — 



2MT 



(1) 



where /(z) = \l{VT+? + z). 

Comparing this background with the number of sterile neu- 
trino events, which is proportional to 0^, we get that the re- 
quired temperature must satisfy 
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>log(l/02) 



Now, if all j3 decay events are considered, the maximal value 
of |q| is Q, and the bound on the temperature is rather rigid 
and reads approximatively, 



T< 



0.7 X 10"^ / m, 
Iog(l/02) Vl keV 



4 /6GeV\ /18.6keVy 
^) \ Q J 



;iK) 



However, one can loosen this bound considerably at the cost 
of the effective source intensity by imposing a kinematic cut 
on the momenta (q)^ < 3MT, which leads to |q|/»Ji ~ in 
^ and to much higher acceptable temperatures 



T< 



log(l/02) VlkeV 



6GeV\ 



M 



)(1K) 



(3) 



Another bound on the experimental sensitivity to the mix- 
ing angle is provided by the requirement that at least one 
decay with sterile neutrino happens during the observation. 
The number of signal events is estimated using the differen- 
tial decay width for massive neutrino dF ^ d^q^E^dq, where 



2me{Q- ^ 



m^) is the electron kinetic energy and 
q— |q| is the neutrino momentum. If the kinematic cut on the 
momenta is much smaller then Q and nis, the number of sterile 
neutrino events can be estimated as 



Vl - ('Ve)A^c, 



(2) 



with A^cut being the number of active neutrino events satisfying 
|q| < C, which is related to the total number of decays A^tot by 

AU _ 35 

wtot"i6Vey 

For the case with generic cut numerical integration of the dif- 
ferential width can be performed. One should also note that 
the cut on the momentum suppresses active neutrino events 
more than sterile neutrino events. 

After this general discussion let us turn to specific num- 
bers. As we see, the best beta decay sources should have small 
mass and small energy release. Then the contribution from the 
thermal motion of the decaying atoms will be small and cre- 
ated sterile neutrino will be not too relativistic, thus making 
it easier to measure its mass. Here the exceptional opportu- 
nity is provided with tritium ^H, which is much lighter than 
all other radioactive elements and has short enough lifetime 
of 12.3 years. However, if the neutrino mass is higher than 
the energy release Q= 18.591 keV in tritium decay other iso- 
topes should be used. The requirements on the temperature 
for the case of ^H2 molecule can be read directly from ( I2I3I ) 
by substituting M by 6 GeV. The bounds (O and A^events > 10 
(3 a with zero background) for tritium are given on the Fig.[T] 
for several different temperatures and exposures. 

One can see, that for T ^ 0.01 K and A^cut ~ 10^^ the exper- 
iment is better than astrophysical X-ray bounds for all achiev- 
able for tritium experiment masses. For a more modest num- 
ber, like A^cut ~ lO'" (i.e. year of observation with 1000 de- 
cay counts per second, coiTesponding to a typical recoil ion 
time of flight in a 1 m long spectrometer) the sensitivity is 
smaller, but still competes with X-ray experiments for a vast 
range of masses. Currently fsT], one can create supersonic 
gas jets with particle densities of about lO" — lO'^ cm^^ and 
temperatures of the order of 0. 1 K and prepare sources with 
magneto-optical traps with densities of lO'" cm^^ and tem- 
perature ^0.1 mK. These methods provide 10^ — 10^ beta 
decays per year for the source size of about 1 mm^. The line 
corresponding to the statistics A^tot ^ 10^ is shown in Fig. [T] 
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together with the hne for 0.1 mK without making any cuts; it 
indicates that the current magneto-optical traps technologies 
for low temperature sources allow to obtain bounds better than 
existing kink searches and to enter in the interesting parameter 
region provided that neutrinos make less than 100% of DM. 

An important point should be emphasized. All the esti- 
mates above have been done with the thermal distribution of 
the source, which is not necessarily the case. As far as the 
number of signal events for keV neutrino is expected to be 
extremely small, the tails of the distribution of thermal veloc- 
ities are important. Large non-Gaussian tails will not spoil 
too much the precision of the mass measurement, but will pe- 
nalize significantly the sensitivity to the mixing angle. This 
may be a problem for some cooling techniques. For example, 
the supersonic jet cooling has some non-Gaussian tails in the 
velocity distribution 1,39^ ,55.1 . 

IV. CONCLUSIONS. 

In this paper we argued that the detailed study of kinematics 
of j3 decays with the help of COLTRIMS may enter into in- 
teresting parameter region for the search of Dark Matter ster- 
ile neutrinos which is complementary to cosmic X-ray mis- 
sions and indispensable for revealing the nature of DM in the 
Universe. Even with currently existing technologies entering 
in the interesting parameter region is possible for light sterile 
neutrinos. Extending the study to compete with X-ray bounds 
for higher masses is a challenging but valuable experimental 
task. 

For a detailed feasibility study of j3 -decay experiments to 
search for DM sterile neutrino a number of extra points, in- 
cluding existence of possible backgrounds, should be clari- 



fied. One obvious background appears from the fact that after 
the j3-decay a fraction (15%) of the (^He^H)+ ions dissoci- 
ates, leading to an error in the determination of momentum of 
the detected recoil ions. However, the momentum after dis- 
sociation is high [44], so that these events do not look like 
the heavy neutrino ones. In addition, one should take into 
account the scattering of ions as a source for possible back- 
ground. Also, a careful analysis may lead to the choice of 
another optimal isotope, which has higher decay energy re- 
lease but is short lived, providing thus larger statistics. A very 
hard problem is the low density of cold atoms (serving as a 
source of beta-decays), available at present. 

A similar experiment can be also made with isotopes de- 
caying by electron capture. In case one has a single-electron 
ion of such an isotope, the final state would contain just two 
particles, and to find the invisible particle mass (v or A^) it is 
sufficient to measure only the energy of recoil ion, which is 
a simpler task than the full momentum reconstruction. How- 
ever, if the initial isotope is not highly ionized, then, gener- 
ally, one or several Auger electrons are emitted after the de- 
cay, carrying away considerable momentum. All these Auger 
electrons should be detected and taken into account (see [4^ 
for ^^Ar case). We leave the comparative study of this type of 
experiments for future discussion. 
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